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Abstract. The problem of linking the structure of a finite linear 
dynamical system with its dynamics is well understood when the 
phase space is a vector space over a finite field. The cycle structure 
of such a system can be described by the elementary divisors of the 
linear function, and the problem of determining whether the system 
is a fixed point system can be answered by computing and factoring 
the system's characteristic polynomial and minimal polynomial. It 
has become clear recently that the study of finite linear dynamical 
systems must be extended to embrace finite rings. The difficulty 
of dealing with an arbitrary finite commutative ring is that it lacks 
of unique factorization. In this paper, an efficient algorithm is 
provided for analyzing the cycle structure of a linear dynamical 
system over a finite commutative ring. In particular, for a given 
commutative ring R such that \R\ = q, where q is a positive integer, 
the algorithm determines whether a given linear system over R n is 
a fixed point system or not in time 0(n 3 log(nlog(g))). 



1. Introduction 

A finite dynamical system is a function / : X — > X, where X is 
a finite set. The dynamics of the system is obtained by iterating the 
function /. Such dynamical systems have a variety of applications, 
such as in engineering, computer science, and computational biology 
[1,3,4]. 

It is a well-known fact in finite field theory that a function / : F™ — ► 
F q , where W q is a finite field of q elements, can be represented by a 
polynomial function. Thus any function / : F™ — > F" can be rep- 
resented by / = (A, . . . , /„), where /< G Fjx x , . . . , x n ). When / is 
a linear system, the dynamics of / can be described using its char- 
acteristic polynomial and minimal polynomial, and the computation 
can be done in polynomial time [1,6]. For general polynomial systems, 
there have been only limited successes in determining the dynamics of 
such systems, except for monomial dynamical systems, where all the 
coordinate functions /j are monomials. 
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In [3], monomial dynamical systems over Z2, i.e. Boolean monomial 
systems, were studied. In [4], the problem of determining whether a 
monomial dynamical system over a finite field ¥ q is a fixed point system 
was reduced to the same question of an associated Boolean monomial 
system and a linear system over a ring of the form Z/(g — 1). In [1], 
the study of fixed point systems was further developed. In particular, 
linear systems were defined for modules over a ring, and a necessary 
and sufficient condition for a linear system to be a fixed point system 
was derived using Fitting's lemma. 

Though the result in [1] does not lead to an efficient algorithm for 
determining whether a linear system over a general finite commuta- 
tive ring is a fixed point system, the computational problem, which is 
ultimately needed in applications, was discussed in some detail in the 
special case where the ring is a finite field, and a computational method 
via the factorization of the characteristic polynomial and the minimal 
polynomial of the linear function was described. As pointed out in [4], 
the approach via characteristic polynomial and minimal polynomial for 
a linear dynamical system over a finite commutative ring faces consid- 
erable difficulties due to the lack of unique factorization (see also the 
comment after Example 4 in [1]). The following example illustrates 
this point. 

Example 1.1. Let f : Z| — > Zg be defined by the 2 x 2 matrix 



Then A k 7^ for 1 < k < 6 and A 6 = 0. Thus f is a fixed point 
system with the only fixed point 0. The characteristic polynomial of A 
is c/iA(A) = A 2 + 6A + 2, which has no root in Z 8 , though A clearly has 
an eigenvector (0,4) T corresponds to the eigenvalue (the eigenvalues 
and eigenvectors of a matrix over a commutative ring are defined as 
usual, see [2]). Note that 



In this paper, we consider a different approach. Our approach is 
based on the fact that there are efficient algorithms for the computation 
of the powers of a matrix: the multiplication of two nxn matrices takes 
at most n 3 operations (the state of the art algorithms use close to n 2 
operations). If A is an n x n matrix, then to compute A m , where m is a 
positive integer, it will take about n 3 log 2 m operations. Therefore, one 
can just work with the matrix of a linear dynamical system directly 
to avoid the difficulties of dealing with the factoring problems over an 
arbitrary commutative ring. 




A 6 = 



(A 2 + 6A + 2)(A 4 + 2A 3 + 2A 2 + 4) (mod 8). 
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In order for this approach to work, one must have a reasonable upper 
bound on the exponents of the powers of the matrix, that is, a reason- 
able upper bound on the number of iterations, that one must compute 
in order to determine the dynamics of a given system. 

Our first observation is, although Fitting's lemma tells us that a 
linear system will be stabilized after a certain number of iterations 
(see [1]), the lemma itself is a fairly general statement: it applies to 
any group G that satisfies both ACC and DCC conditions on normal 
subgroups and any normal endomorphism / of G (see [7, p. 84]). While 
for the systems that we are interested in, the groups involved are finite 
abelian groups, and therefore, we should be able to derive more precise 
information on how many iterations it will need in order for a given 
system to reach a certain type of stabilization status. 

Our second observation is, the upper bound on the number of itera- 
tions also depends on the size of space. This can be seen from Example 
1.1, where it takes 6 iterations for the system to be stabilized. This 
can also be seen by just considering the simplest type of linear sys- 
tems on Tig, where q is a positive integer, namely the ones defined by a 
scalar multiplication. For such a system, the matrix size is 1, but the 
dynamics of the system depend on q. If the system is defined by the 
multiplication of an element 1 < a < q, then one either needs to know 
the prime factorizations of a and q or needs to compute the powers of 
a to derive the dynamics of the system. Therefore, certain assumption 
on the size of q must be made. Here we assume that the size of q is 
comparable to the size of any integer that we maybe able to factor in 
the foreseeable future. We believe it is reasonable to make this assump- 
tion. Under this assumption, the numbers log 2 q and log 2 (log 2 g) are 
relatively small: the RSA keys are typically 1024 — 2048 bits long and 
log 2 (log 2 2 2048 ) = ll. 

This paper is organized as follows. In Section 2, we develop the basic 
theory that lays the foundation for an efficient algorithm. In section 3, 
we describe an algorithm for determining whether a linear dynamical 
system over a finite ring is a fixed point system or not and give two 
examples of linear fixed point systems over finite rings which are not 
fields. In Section 4, we conclude with some discussions and an example. 
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2. Main results 

Let R be a finite commutative ring with q > 1 elements. Let the 
prime factorization of q be 

t 

i=i 

We shall view the elements of R n , where n is positive integer, as 
column vectors, and denote by e^, 1 < i < n, the canonical basis (if R 
has 1). For a function / from a set to the same set, we write 

/m = /o/o ... o/ 

V «, ' 

to copies 

if to is a positive integer. If to is a positive number, not necessary an 
integer, then by writing f m we mean f\ m \ where [to] is the smallest 
integer greater than or equal to to. Our first theorem upper bounds 
the number of iterations needed for a linear system to reach a certain 
stable status. 

Theorem 2.1. Let n be a positive integer, and let f : R n — > R n be a 
linear function. Then for any nonnegative integer k, we have 

jn\og 2 (q)+k _ ynlog 2 (?)^n^ 

If R is a field, then the factor log 2 (g) is not needed, that is 

f n+k (R n ) = f n (R n ). 

Proof. We first consider the general case when R is a commutative 
ring. View R n as an /-module, set M = R n , and consider a sequence 
of /-submodules of M defined by 

(2.1) M DM 1 = /(M ) D ■ ■ ■ D M r = f(M ) D ■ ■ ■ . 

Since each M r (r > 0) is a finite abelian group and 

t 

\M \ = q n = Hpf% 
i=i 

by Lagrange's theorem, we have 

t 

|M r |=IJp?, 

i=\ 

where < r\ < nti. Thus, if M r ^ M r+1 , then 

|M r+1 | < \M r \/ Pi 
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for some 1 < i < t. Therefore either there is an 

t t 
(2.2) r < nti = n ^ tj := s, 

i=i i=i 

such that M r = M r+1 , or we must have \M S \ = 1. In any case, f(M s ) = 
M s . Since s < nlog 2 (g), the first statement follows. 

If R is a field, then the modules M; are vector spaces over R. So 
if Mi D M m , then dimM i+ i < dimM; - 1. Since dimM = n, the 
desired result follows. □ 

Next, we give a general lemma about fixed point systems on a finite 
set. We remark that one can almost read out the proof of the lemma 
from the proof of Theorem 2 in [1]. Here we give a proof which sheds 
some light from a different view. 

Lemma 2.1. If X is a finite set and f : X — > X is a function such 
that f(X) = X , then f is a fixed point system if and only if f is the 
identity function. 

Proof. Since X is a finite set, f{X) = X implies that / is also injective. 
Thus / is a permutation of the set X. Writing / as a disjoint product 
of cycles, we see immediately that / is a fixed point system if and only 
if all the cycles have length one, that is, / is the identity function. □ 

Recall that an element a in a commutative ring R with 1 is called 
a unit if it is invertible. The following is an immediate consequence of 
Lemma 2.1. 

Corollary 2.1. Let R be a finite commutative ring with 1. Let A : 
R n — ► R n , where A is an n x n matrix over R, be a linear dynamical 
system. If A ^ I and det A is a unit in R, then A is not a fixed point 
system. 

Now we give a criterion for a linear dynamical system over a finite 
ring to be a fixed point system. 

Theorem 2.2. Let R be a finite commutative ring with q elements, let 
n be a positive integer, let f : R n — > R n be a linear system, and let 
A be the matrix of f with respect to the canonical basis (if R has 1). 
Then f is a fixed point system if and only if f nXo ^ii)+ 1 — f nl °S2(i) ^ or 
equivalently A nlo ^ +1 = A nXo ^ q \ If R is a field, then the condition 
simplifies to f n+1 = f n or A n+1 = A 11 . 

Proof. By Theorem 2.1, 

j^jnlog 2 (g)^n^ _ jn\og 2 (q) 
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Thus, by Lemma 2.1, / is a fixed point system if and only if 
which is equivalent to 

That is / rilo §2(9)+ 1 — f»Hog 2 (g)_ □ 

Theorem 2.2 provides an efficient algorithm to determine whether a 
linear dynamical system over a finite commutative ring is a fixed point 
system, which will be discussed in the next section. The results in this 
section also reduce the study of a general linear dynamical system over 
a finite commutative ring to an invertible non-fixed point system. 

3. Algorithms and Examples 

In this section, we first describe an algorithm based on Theorem 2.2 
for determining whether a linear system A : Z™ — > Z™ is a fixed point 
system or not, where q > 1 is an integer and A is taken to be the form 
of an n x n matrix. We choose Z g as the base ring for the simplicity of 
the statements, the same algorithm works for any ring of the type 

Z ?1 X Tj q2 

x ••• x Z, fc , 

as well as for any finite commutative ring with 1 as long as the opera- 
tions of the ring are implemented. 

The algorithm is called an LFPS (Linear Fixed Point System) test. 

Algorithm 3.1. LFPS test. 

Input: Two positive integers n and q > 1, an n x n matrix A over 
Z q , and bt-il 1 ^ 1 + 6 t _ 2 2 < ~ 2 ■ ■ • + b\2 + bo, the binary representation of 
\n\og 2 q\. 

Output: true or false. 

(1) X^I 

(2) for % from t — 1 down to do 

X ^XX 
if bi — 1 then 
X <— AX 

(3) if X = XA then 

return true 
else 

return false 
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Let us explain this algorithm in more detail. In step (1) the (ma- 
trix) variable X is initialized by the identity matrix /. The main com- 
putation of A^ nlog2q ~\ is performed in step (2) using the "square and 
multiply" method. Since 6 t _i = 1 (the leading bit of |~nlog 2 g~|), at 
the beginning (i.e., % — t — 1), X first becomes II = I, then becomes 
X = AI = A. After this, for each i with t — 2 > i > 0, the value in 
X becomes the square of the value previously stored in X. If hi = 1, 
then the value of X is further updated to be the product of A and the 
previous value. At the end of step (2), the value in X is A^ nlog2q \ For 
example, if A is a 6 x 6 matrix over Z 3 . 7 , then |~61og 2 21] = 27 and by 
the "square and multiply" method: 

A 27 = ^1.24+1.23+0.23+1.2+1 = ^ A ) 2 A) 2 JaJa. 

In step (3), the result of Theorem 2.2 is applied. Since the value of X 
is now A^ nlog2q \ the system is a fixed point system if X = XA, and 
the program returns true; otherwise, the system is not a fixed point 
system and the program returns false. 

Suppose two matrices over Z g can be multiplied with 0{n u) ) op- 
erations, by using Strassen's algorithm, u < log 2 7. This number 
can be further reduced, see [5]. The cost of running LFPS test is 
0(n UJ (\og 2 n + log 2 log 2 g)). Under our assumption that log 2 log 2 g is 
small, determine whether a linear system over a finite ring is a fixed 
point system or not can be done with 0(n 3 ) operations. If R is a field, 
then the number of operations required is 0(n UJ \og 2 n). 

As long as the problem of determining whether a linear system is a 
fixed point system is concerned, a comparison of the computational cost 
analysis given in [1] with the analysis given above shows, in addition to 
its simplicity, that our algorithm is at least as efficient as the approach 
via the characteristic polynomial and the minimal polynomial even for 
the case of finite fields. 

Next we give two examples of fixed point linear systems over finite 
rings. The first example is over the ring Z 2 4. 

Example 3.1. The system A : Z 24 — > Z 24 defined by 

/15 7 7 l\ 

7 117 

7 7 7 11 ' 
\U 8 15 6 J 

is a fixed point system. This can be verified by using Algorithm 3.1 to 
compute A 41oS2 24 = A 16 (A iterations) and verify that A 16 = A 17 . The 
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"stabilized" matrix is 



16 



(12 


1 


2 


ll\ 





4 


8 


12 


4 


3 


6 


1 


V 12 


1 


2 





A 



The second example describes a fixed point system over the ring 
Z 3 2. 5 . 

Example 3.2. T/je system A : Zg 2 



A : Z^ 2 . 5 




Zg 2 5 defined by 


/36 23 


32 


9\ 




27 32 


30 


25 




32 25 


13 


28 




\32 8 


41 


40 y 





/O 


9 


9 


27\ 


10 


27 


12 


26 


35 


18 


33 


19 


V 5 


27 


42 


3l) 



is a fixed point system. The "stabilized" matrix is 



^4riog 2 (3 2 -5)l = A 24 = 



We remark that the number r such that A r = A r+1 can be smaller 
than our theoretical bound n\og 2 q in some cases. In the second ex- 
ample above, r = 6, i.e., we have A 6 = A 7 . Our algorithm can be 
refined so it terminates before the iteration process reaches the theo- 
retical bound if r is small enough, say r < a/ n log 2 q. But we believe 
that the gain is not significant by doing so. 



4. Conclusions 

We have provided an efficient algorithm to determine whether a lin- 
ear dynamical system over a finite commutative ring is a fixed point 
system. As an application, our result together with the results in [3] 
and [4] should settle the problem of determining whether a monomial 
dynamical system over a finite field is a fixed point system. 

When the system is not a fixed point system, a natural problem is 
finding the cycles of the system. If R is a field, then under the assump- 
tion that the elementary divisors of a linear system and their orders 
(the order of a polynomial g is the least positive integer k such that g 
divides x k — 1) can be computed efficiently, the cycles can be computed 
by a theorem due to Elspas (see [6]). Obviously, the implementation 
via such approach is quite involved, in particular the computation of 
the orders of the elementary divisors. The orders of the elementary 
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divisors are the lengths of the cycles. If the lengths of the cycles can 
be found, then the cycles can be obtained. For example, suppose that 
/ is linear dynamical system over a finite commutative ring R with 1, 
and suppose that the lengths of its cycles, say 

= k < ki < . . . < k m , 

are known. Then the cycles can be computed by solving the linear 

systems: 

(A k * -I)X = 0, 0<i<m. 
Example 4.1. Consider the system A : Z\ 05 — ► Zf 05 defined by 

[70 27 5 26 \ 

35 98 104 99 

81 85 78 102 ' 

\27 97 13 69 / 

Since detA = 2 (mod 105) is a unit in Z 105 , Corollary 2.1 implies 
that A is not a fixed point system. Since A 24 = I and A k ^ I for 
< k < 24, we see that the cycles lengths are the factors of 24. With 
some computation, one can find the cycles lengths, they are 1,2,24. 
The only cycle of length 1 is 0, there are 5512 cycles of length 2, and 
5064150 cycles of length 24. 

However, the search for the cycle lengths seems to be exponential. 

Computations of linear systems over finite commutative rings are 
basic, since one typically handles the other computational problems by 
reducing them to the ones about linear systems, and for systems over 
finite fields, the reduction can result in linear systems over commutative 
rings which are not necessary fields. Developing efficient algorithms 
over commutative rings deserves further attention (see also [1]). 
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